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ABSTRACT 

Two  kinds  of  finite-element  analyses  are 
presented  for  the  investigation  of  the  large 
deformation  phenomenon  in  laminated  composite 
structures,  especially  shells.  The  first  kind  of 
finite-element  analysis  employs  the  general 
Incremental  variational  formulation  as  well  as  the 
total  Lagranglan  description  of  motion*  The  element 
adopted  Is  the  “degenerate"  three-dimensional 
element*  The  secoad  analysis  adopts  a  formulation 
based  on  deformable  shell  theory,  and  the  plate- 
bending  element  is  used.  Numerical  results  for 
bending  are  presented  for  plate  and  shell  structures 
of  isotropic  as  well  as  orthotropic  co^oeltlon.  The 
resuts  are  in  good  agreement  with  results  available  in 
the  literature. 

INTRODUCTION 

The  shell  is  a  common  load-carrying  element  in 
most  Industrial  equipment,  especially  in  aerospace, 
nuclear  and  offshore  engineering.  Although  linear 
models  suffice  In  many  aspects  of  engineering,  a 
linear  analysis  of  shells  yields  results  too  gross  and 
Inaccurate  to  be  useful.  The  nonlinear  behavior  of 
structures  should  be  taken  Into  account  If  e  thin 
shell  structure  design,  dlctatad  by  economic 
considerations,  is  desired.  There  exist  a  number  of 
shell  finite  elements  in  tha  literature.  A  discussion 
of  various  alemants  is  glvsn  by  Gallsghar  [1]. 

Tha  degenerate  three-dlmaas ional  element, 
originally  proposed  by  Ahmad,  Irons  and  Zlenklevlcs 
[2]  for  linear  analysis  of  thin  and  thick  shells,  was 
applied  to  problems  dealing  with  the  geometrical  and 
material  nonlinearl t lee  la  [3]  ead  [4].  It  la  now 
believed  that  this  element  performs  well  In  shell 
analysis.  Worssk  [3]  proposed  e  bilinear  degenerate 
element,  ead  Worsak,  at  el.  [6]  developed  an  eight- 
node  element  that  has  relative  displacement  degrees - 
of-freedom  In  place  of  the  cumbersome  finite 
rotations.  Tha  derivation  of  both  elements  was  based 
on  the  concept  of  the  degenerate  3-dlmenalonal  element 
model.  Gallagher  (7),  end  Brabble  ead  Connor  (8] 
extended  this  concept  to  the  geometrically  nonlinear 
analysis  of  shells.  The  degenerate  3-D  element  la 
usually  accompanied  by  incremental  form  let  lone  in  the 
analysis  of  geometric  end  materiel  nonlinear  it  lee  In 
problems  of  structures ;  various  lacremsatal 
variational  fonmtletlons  for  the  total  Lagranglan  as 
well  as  tha  updated  Lagranglan  approach  are  presented 
in  (f,10J. 

In  modeling  laminated  composite-material  shells, 
Domg,  Hater  mad  Taylor  (11)  formulated  e  theory  of 
thin  shells  laminated  of  anisotropic  notorial,  which 
la  as  extension  of  the  theory  developed  by  Stavefcy 
112)  for  laminated  anisotropic  plates  baaed  oe  Love’s 
first  approximation  theory  of  shells,  hrafcateeh  end 


Rao  [13]  used  a  doubly  curved  quadrilateral  element  to 
model  laminated  shells  of  revolution;  there  are  12 
unknown  displacements  at  each  node.  Noor  (14-16]  used 
a  mixed  formulation  to  compare  several  alemants  in  tha 
linear  analysis  of  laminated  shells.  Hsu,  Reddy  end 
Bert  (17]  employed  the  sheer-deformable  element 
proposed  by  teddy  (18]  to  analyte  shells  composed  of 
bimoduler  materials.  Noor  (19)  employed  the  classical 
shall  element  that  Includes  the  geometric  nonlinearity 
to  solve  e  spherical  shell  problem.  A  degenerate 
element  associated  with  the  updated  Lagranglan 
formilation  la  presented  In  (20]  to  solve  e  cylinder 
problem. 

The  present  paper  teals  with  the  geometrically 
nonlinear  analysis  of  layered  composite  shells,  using 
e  3-D  degenerate  element.  For  the  sheer  deformable 
shell  element  description,  see  Reddy  [21,22]. 

A  3-D  DEGENERATE  KLIMENT  FOR  LAYERED  COMPOSITE  SHELLS 

Consider  the  notion  of  e  material  body  In  a 
Euclidean  apace,  aa  shown  in  Figure  1.  Here  Cc 
denotes  the  configuration  of  the  body  at  time  t.  To 
facilitate  the  description  of  the  kinematics  of  tha 
continuum  at  any  time,  the  initial  state  of  the  body 
can  be  chosen  ee  the  reference  configuration  to  which 
ell  subsequent  configurations  are  referred  to;  this  Is 
known  ee  the  total  Lagranglan  description.  If  tha 
coordinates  are  updated  at  each  Incremental  time-step, 
and  tha  kinematics  of  the  continuum  la  described  in 
tarns  of  the  current  configuration,  tha  description  is 
celled  the  updated  Lagranglan  description.  It  is 
noted  that  e  load  step  is  used  instead  of  the  time 
step  If  we  ere  solving  the  static  bending  problem. 


Figure  1.  Description  of  the  notion  of  e 
a  continuum 


Invoking  the  principle  of  virtual  displacements, 
one  enn  express  the  equilibrium  equation  of  the  body 
et  time  t  4-  At  me: 


t4dt 

o 


(i) 


t4At 


6V 


where  the  left  super script  refers  to  the  configuration 
st  which  the  quantity  is  calculated;  the  left 
subscript  is  the  configuration  to  which  the  calculation 
la  referred,  la  the  second  Piola-Ki rchbof f  stress 
tensor  and  c. .  is  the  Green-Lagranglaa  strain  tensor 
defined  as  J 


O  1J  V  o  t.J  ouj,l 


and  denote  the  component*  of  surface  and  body 
forces,  and  V0  and  Ag  the  volume  and  area, 
respectively,  of  Che  body  in  the  Initial 
configuration.  Using  the  relations 


oSlj 


«S13 


+  oSiJ 


(3.) 


and 


t+At 


o"l  +  oV 


(3b) 


where  and  Qui  are  the  Increaental  stress  and 
displacement  components,  we  write 


t*tc 

•  ij 


+  oeij  * 


where 


(4) 


0*1J  “  o*iJ+  onij 

0*1J  ”  'J^oui,J+ouJ,i4ouk,l  ouk,J+oak,i  o“k,J ^ 
ontJ  “  I  o°lt,i  o“k,j 

The  linear  constitutive  relatione  are  given  by 
o*lj  ■  Cijk»  o4kl* 

The  equilibrium  equation  (1)  then  becomes 


(3a) 

(3b) 

(3c) 


(6) 


CtK^l  ♦  (Km.))('»l  *  -  1*1  ♦  1*1  (•) 

vh«r«  (K^J ,  [K^J  ud  |k)  ir<  th*  Uuar  stiff.*., 
matrix,  nonlinear  atlffneaa  eatrlx  and  unbalanced 
terms,  rsspsctlvely;  {?}  Is  ths  forcing  term  from 
external  loading.  To  solve  these  equations,  ths 
Mswton-isphson  or  ths  modified  Wswton-taphaon  method 
can  be  used.  If  the  ]  le  updated  after  much 
iteration,  one  has  the  newton  -laphaoo  method.  When 
the  [Kw]  is  kept  constant  during  each  step,  one  has 
the  Modified  Wswton-taphaon  method. 

mnSKICAL  XK  SUITS 


The  f ini ts-elemsn t  formulation  preaantad  ha rein 
Includes  the  treneverqe  sheer  effect,  although  we  have 
assumed  that  the  straight  line  normal  to  the  mid-plena 
before  deformation  remains  straight  after 
deformation.  The  reduced  Integration  technique  la 
uaed  In  tha  numerical  integration,  bie  to  tha  biaxial 
symmetry  In  the  problems  considered »  only  e  quarter 
domain  le  considered.  A  2x2  meeh  of  8-node  end  9-node 
elements  era  employed  In  this  analysis. 


PKOBLEM  1:  Laminated  Composite  Plates 

A  simply  supported,  2-layer  crossply  square  plate 
made  of  graphite-epoxy  materiel  subjected  to  a 
uniformly  distributed  load  le  analysed.  The  material 
properties  ere 

Kj/Ej  -  40  ,  C12/K2  -  0.5  »  vJ2  -  0.25  (9) 

end  the  boundary  conditions  ere 


w  ■  ♦  •  0  st  y  -  0,a, 
w  •  ty.  -  0  et  x  -  0»a, 


(10) 


The  numerical  solutions  for  both  elements  ere 
presented  lu  Table  1.  It  le  deer  that  the  DST 
predicts  lower  deflections  then  the  3-D  theory. 


Table  1.  Comparison  of  the  nonlinear  center  - 

deflection  (w/h)  of  Problem  1  obtained  by 
the  3-D  degenerate  end  shell  elements. 


P  1.0  2.0  3.0  4.0  5.0  6.0 


3-D  1.0414  1.6431  2.0604  2.3871  2.6595  2.8955 
DST  1.0754  1.6042  1.9971  2.3008  2.5560  2.7716 

P  -  (P0(./h)V«2J  s  10-2.  ~ 

3-D  •  degenerate  3-D  element 

DST  •  deformable  shell  (2-D)  theory  (21] 


W. 


.CiJk*  o*ki  o**tj  iyo  +1  oSiJ  o4* 


dV 


t44t. 
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(7) 


which  is  s  nonlinear  equation  la  terms  of  (he 
incremental  displacements  u,  v  and  w.  For  the 
degenerate  element  (see  (23])  this  equation  cun  be 
cast  In  tens  of  the  middle  surface  nodal  variable# 
*o»  *•’  V  •  sad  6.  In  matrix  form,  me  have 


PtOSLSM  2s 


Isotropic  Cylindrical  Penal  With  Uniform 
Loading 


Consider  the  circular  cylindrical  panel  shewn  la 
Flgmra  2.  Tha  paaal  le  clamped  along  ell  four  edges, 
tad  le  subjected  to  uniform  normal  loading.  A  3x3 
aaah  of  9-ooda  elements  Is  used.  Ths  dlmtaaioaal 
csnter  deflections  obtained  by  both  elemeate  ere 
listed  In  table  2,  Once  again  we  note  that  the  DST 
element  predicts  lower  deflections. 


•  Codes 

*  'n*.  or 
?'  c  J  o  1 


^"L  *■*!•— i 


R*2S40  mm 
Ls254  rrm 
h*6.35  mm 
0*0. 1  rad# 

E»3. 10275  m/m2 
\>*0.3 


Figure  2.  Cylindrical  panel  used  in  Problems 
2  and  3 


Table  2.  Comparison  of  Cha  nonlinear  canter  - 

deflection  w  (in  am)  of  Problem  2  obtained 
by  the  3-D  degenerate  element  and  shell 
element 

P  "  1*0  lTs  1/75  £5  2725  2.50 

(kM/m2) _ 


DST  1*3827  3.8828  5.9449  7.5132  8.546  9.2949 


h/L  -  0.01 

Figure  3.  Spherical  shell  used  In  Problem  4 


Table  4.  Comparison  of  the  nonlinear  center - 

deflection  (v/h)  for  Problem  4  obtained  by 
various  Investigators 


p 

2  4 

5 

6  8 

10 

DST 

0.70918  1.8591 

_ 

2.8299  3.6260 

3.8804 

3-D 

0.69053  1.7607 

- 

2.7649  3.4170  3.8897 

MOOR  0.6806  1.741 

2.305 

3.436 

3.916 

PROBLEM  5:  Laminated  Cylindrical  Panel  VI th  Uniform 


PROBLEM  3;  Hinged  Circular  Cylindrical  panel  With  Pin 
Load 

Using  the  same  geometry  of  the  panel  in  Problem 
2,  but  with  thickness  12.7mm  and  boundary  conditions 
such  that  the  curved  edges  are  hinged  and  Che  straight 
edges  are  free,  A  concentrated  normal  loading  la 
applied  at  the  center.  Table  3  contains  the  numerical 
results: 

Table  3:  Comparison  of  the  nonlinear  center- 

deflection  w  (in  mm)  in  Problem  3  obtained 
by  the  3-D  degenerate  element  and  shell 
element 


0.5  1.0 


2.0  2.2 


3-D  1.3784  2.9786  4.9390  5.8965  7.719  9.7715 
DST  1.3707  2.9591  4.9088  -  7.627  9.4931 


PROBLEM  4 :  Laminated 
Uniform  Lo 


oslte  Spherical  Shall  With 


Cons Ida r  a  9-layered  cross-ply  laminated  shell  of 
graphite-epoxy  (see  Figure  3).  The  four  edges  are 
simply  supported;  i.e., 

u  *  t  ■  ♦  *  0  it  y  •  0,  a, 

ui> 

f  •  »  ■  9  •  0  it  x  -  0,  a. 

A  uniform  inward  loading  is  applied  on  the  shall.  Tha 
numerical  results  are  tabulated  In  Table  4.  The 
present  results  compere  very  sell  with  those  of  Moor 
(19).  In  this  case,  the  DST  element  predicts  larger 
deflections  while  Moor’s  mixed  element  predicts  lower 
deflections  than  thosa  predlctad  by  the  3-D  degenerate 
element. 


The  same  shall  gaomstry  aa  In  Problam  3,  but  with 
2  layers  f0*/90*]  Is  used.  The  material  is  graphlte- 
apoxy  and  all  four  edgaa  ara  simply  supported,  i.e., 

u  •  »  •  *  •  0  it  y  ■  0,a, 

*  (12) 
0,a. 

The  results  ere  compared  In  Table  3. 

Table  5.  Comparison  of  the  nonlinear  center - 

def lection  (w/h)  of  Problem  5  obtained  by 
the  3-D  degenerate  element  and  shall 
element • 


-  V 

P*  (— tt)  50 

V 


100  200  300  350 


DST  0.6326  1.2533  1.9739  2.3988  - 
3-D  0.6014  1.2340  1.9812  2.4317  2.471 
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ABST1ACT 

Results  of  a  preliminary  study  of  the  three- 
dimensional  ,  geometrically  nonlinear,  finite-element 
analysis  of  tha  banding  of  laminated  anisotropic 
composite  platan  ara  presented.  Individual  lamina « 
ara  assumed  to  ba  homogeneous,  transversely  Isotropic, 
and  llnaarly  elastic.  Tha  investigation  uclllsas  a 
fully  three-dimensional  Isoparametric  finite  element 
with  sight  nodes  (l.e.,  linear  element)  and  24  dagraas 
of  freedom  (thraa  displacement  coaponants  par  node)* 
Numerical  rasults  of  tha  linear  analysis  ara  compared 
with  tha  exact  solutions  obtained  by  Psgeno.  Tha 
present  rasults  converged  to  tha  exact  solution  as 
aesh  is  refined.  Due  to  lack  of  space,  rasults  of  tha 
nonlinear  analysis  ware  not  included  hare;  tha  results 
will  appear  elsewhere. 

introduction 

Tbs  development  of  slapls  ft nits  slsasnts  for  tha 
analysis  of  plates  has  bean  a  subject  of  continued 
Interest  among  researchers  of  computational  solid 
mechanics.  Most  of  tha  currently  available  slsasnts 
ars  bsssd  on  two-diatnslonal  thln-plate  theory  In 
which  tha  plana  sections  normal  to  tha  mldsurface 
before  deformation  ara  assumed  to  remain  straight  and 
normal  to  tha  mldsurface  after  deformation  (!•«.» 
transverse  shear  strains  ara  assumed  to  bt  aero)  sad 
thlck-plst*  theory  that  accounts  for  constant  or 
linear  variations  of  trensverse  shear  strains  through 
the  thickness.  With  the  Increased  use  of  advanced 
fiber  composite  materials  in  aerospace  sod  automotive 
structures,  s  mors  accurate  prediction  of  stresses  in 
composite  plates  is  needed.  Composite  materials 
exhibit  higher  st if fosse -to- weight  ratios  and 
lncraaaod  corrosion  resistance  compared  to  isotropic 
materials.  The  anisotropic  material  proport las  of 
layered  compos 1 too  can  bo  varied  by  varying  tha  fibre 
orientation  and  stacking  sequence.  Wills  this  feature 
gives  tha  designer  an  added  degree  of  flexibility,  the 
stiffness  mismatch  of  tha  orthotropic  layers  bonded 
together  with  different  fiber  orientations  loads  to 
interlaminar  stresses  in  the  vicinity  of  free  edges. 
For  certain  stacking  sequences ,  loading,  and  boundary 
conditions  Interlaminar  stresses  con  bo  so  largo  that 
they  dictate  the  design  of  the  structure.  An  accurate 
prediction  of  the  interlaminar  stresses  Is  possible 
only  when  s  three-dimensional  theory  is  used. 

Previous  analyses  of  layered  composite  plates  can 
ba  divided  Into  two  groups:  (l)  analyses  based  on  o 
lamination  theory,  and  (11)  analyses  based  on  a  fully 
three-dimensional  theory.  Tha  lamination  theory  is  an 
extension  of  the  classical  plats  theory  (CFT)  or  the 
Releener-Mtfidlln  shear-deformable  plats  theory  (SOT) 
to  layered  composite  plates.  The  first  leal  nation 
theory  is  apparently  due  to  teleener  and  Stsveky  (1), 
Tang,  Morris,  end  Stsveky  [2]  presented  a 
generalisation  of  the  Selaener-Ktsdlln  chick-plate 


theory  for  homogeneous ,  isotropic  platss  to 
arbitrarily  lamina  tad  anisotropic  plstss.  Whitney 
end  Pagano  (3]  (also  sss  Reddy  and  Chao  14] )  presented 
closed-form  solutions  to  the  theory  whan  applied  to 
certain  cross -ply  and  angla-ply  rsctangular  platss. 
lad dy  [5]  presented  a  finite-element  analysis  of  tha 
lamination  theory.  A  higher-order  lamination  theory 
that  accounts  for  a  cubic  variation  (as  opposed  to 
linear  in  (2-5])  of  the  lap lane  displacements  and 
quadratic  variation  of  the  transverse  displacement 
through  tha  thickness  is  presented  by  Lo,  Christensen, 
and  tfo  (61,  and  hybrid-stress  finite -element  analysis 
of  tha  theory  la  presented  by  Spllker  (7].  A 
generalisation  of  tha  von  Kerman  nonlinear  plate 
theory  of  isotropic  platss  to  anisotropic  plstss  is 
due  to  Ebcloglu  (8].  Chla  and  his  colleagues  (9,10) 
presented  approximate  numerical  solutions,  using  the 
perturbation  method  and  tha  Galerkln  method,  of  tha 
von  Kerman  aquations  associated  with  layered  composite 
plates.  Recently,  Reddy  and  Chao  (11,12]  presented 
finite-element  analyses  of  the  same  aquations. 

In  the  lamination  theory  it  is  assumed  that  the 
laminate  la  In  a  atata  of  plana  a trees  (an  assumption 
carried  from  tha  plate  theory)  and  integrals  through 
the  thickness  of  a  laminate  are  equal  to  tha  sum  of 
integrals  through  the  thickness  of  individual 
laminae.  These  assumptions  lead  to  inaccurate 
prediction  of  Interlaminar  straaaas  at  the  free  edges; 
away  from  free  edges  the  lamination  theory  solutions 
ars  vary  accurate.  Thus  the  lamination  theory  la  not 
accurata  In  a  boundary  layer  region  which  extends 
Inward  from  tha  edge  to  a  distance  approximately  equal 
to  tha  laminate  thickness • 

The  fully  three-dimensional  theory  la  nothing  but 
tha  linear  elasticity  theory  of  a  three-dimensional 
solid.  Ft pas  and  Pagano  (13]  and  Pipes  (14)  used  a 
finite-difference  technique  to  solve  the  quael-three- 
dlnenslonal  elasticity  aquations  for  laminates  (also 
sea  Hsu  and  Varakovlch  (XS1).  Us  [161,  Dana  [17], 
and  Dana  and  Barker  [18]  used  a  cubic,  three- 
dimensional,  Isoparametric  element  with  72  degrees  of 
freedom  to  analyse  laminated  plates.  Tha  numerical 
results  in  these  studies  agree  very  well  with  those  of 
Pagano  (19,20/.  The  studies  were  confined  to 
geometrically  (aa  wall  aa  materially)  linear 
analyse*.  Recently,  Griffin,  Kimat,  sad  Bsrakovleh 
[21]  Investigated  three-dimensional  inelastic  (l.o., 
notorial  nonlinearity)  finite-element  aaalyoia  of 
laminated  composites,  while  Brockaen  [22]  developed  a 
finite-element  program  (called  MAGMA)  for  three- 
dimensional,  geometrically,  nonlinear ,  analysis  of 
Isotropic  plates. 

The  present  study  Is  motivated  by  tha  lack  of 
finite-element  results  for  three-dimensional , 
geometrically  nonlinear,  saalysls  of  layered  composite 
plate* •  The  present  paper  gives  e  finite -element 
fomelatlen  of  tbs  nonlinear  aquations  governing  a 


Hourly  elastic,  3-D  continuum.  Numerical  results  of 
the  Ilnur  analysis  nr*  presented  for  several  problur 
and  the  lluar  formulation  la  validated  by  comparing 
th*  prasant  raaulta  with  thou  of  Pagano  [19,20]  and 
others.  Raaulta  of  th*  nonlinear  analyala  will  appear 
elsewhere. 


Cooaldar  a  laminate  <fl)  composed  of  If  orthotropic 
layaraa  with  axes  of  alaatlc  symmetry  parallel  to  the 
Plata  axaa.  The  laminate  la  auhj acted  to  normal 

traction  t^  •  qlx^xj)  at  lta  upper  surface  (t.e., 

•  h/2).  The  constitutive  aquations  for  any  layer  are 
given  by 


**11  **12  **13  0 

**12  **22  **23  0 

**13  **23  **33  0 


**»*  0 


given  In  terms  of  via  Rqe.  (1)  and  (3): 

°1J  “  7 \j*ilu«,o+Un,«+Uk,M“k,a1  <4) 


**1J  .***-«.  J  -  ■ 


ijmn  0  ,  otharvlaa. 


Nut,  we  assume  that  the  displacements  uB  are 
Interpolated  by  expressions  of  the  form. 


um  "  1  Vo  {m  *  X*2*3) 

o«l 


where  #a(x,y,a)  (a  -  1,2»«..,S)  ere  the  trlllnur 

Interpolation  functions  of  the  eight-node 
Isoparametric  a lament  in  thru  dimensions. 
Substituting  tq.  (6)  Into  (5),  we  obtain 


1  1  c*-1*3*3;  . *)»  <*) 

nel  a«l 


and  the  governing  field  equations  are  given  by  (In  the 
absence  of  body  forces) 


where  o^  and  are  the  rectangular  Cartesian 

components  of  the  stress  and  strain  tensors, 
respectively,  and  are  the  material  properties  of  s 

lamina  in  the  laminate  coordinates.  In  (2)  the 
summation  convention  on  repeated  Indices  is  used.  The 
strain-displacement  equations  are  given  by 


cu  *i  l*i,j  +“j,i  ♦  ViVj1 


To  complete  the  description  of  the  field 
equations,  Eqs.  (l)-(3)  should  be  adjoined  by  boundary 
conditions.  At  any  point  of  the  boundary  of  the  plate 
one  should  specify  ou  of  the  following  two  types  of 
boundary  conditions : 

(1)  essential  (or  geometric)  boundary  conditions 


(11)  natural  (or  dynamic)  boundary  conditions 

*.*  Vu<#i  +  Vi)*  V  <*>b 

Finite-Element  Formulation 

Here  we  present  s  displacement  finite-element 
model  of  the  equations  (l)-(A).  To  this  sad  we 
construct  the  variational  formulation  of  the  equations 

over  an  arbitrary  element  of  the  finite-element 
uah.  We  have  (see  (23,  p.  362]) 

**  ,  J4"«.J,V4*l'Vi),*,l4*24*3  *  1  ,  '.4V 


W  („  V-  *  £  <  ^>^lte2ta3 


£  -  /  t  *  d*  .  <J,«,a  -  1,2,3)  (10) 

*  rC)  *  * 

E»«ry  l«opar*Mtrlc  ft  ale*  ilnnt  of  th* 

flnlte-elemrat  suh  can  be  geureted  from  the  muter 
element  (which  le  e  cube  of  side  2)  via  the 
transformation  (am  figure  1) 

*1  "  *l^lf^2f^3^ 


-  t  *iV€H€2t€3)  (1  "  lb2’3)  <U) 

where  x*  are  the  global  coordinates  of  the  element 

nodes.  Therefore,  the  Integrals  In  Iq.  (10)  cu  be 
transformed  to  the  muter  e lament  and  evaluated 
numerically  using  the  Gauss  quadrature.  The 

transformation  of  to  **—  Is  performed  u  follow# 
<su  teddy  (24])i  1  *1 


-  ur* ^ 


-  <d.t  m)«t«2«3 


lm  <  dmotm  th*  T*rl*tloa*l  tyabol,  and  t* 
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*2 

*»2 
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rq  — 
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For  example,  consider 

»♦.  »♦« 
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i  i  i 


E  1  £  ^  (PT*pj*Wjwr 
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(13) 


where  and  Wj  trt  the  Gauss  point*  and  weights, 
respectively,  and  tba  In  tag  rand  F*jj  (i,J  • 

1,2,3;  o,6  -  1,2, •••,8)  la  given  by 


Jkl  balng  the  elements  of  tha  Invar an  of  tha  Jacobian 

matrix,  [ J] .  Thla  procadura  can  ba  Implanantad  on  a 
digital  computer,  and  tha  alamant  coefficient  matrlcaa 
In  Eq«  (10)  can  ba  evaluated  numerically*  Tha  lxlxl 
Ganna  quadrature  was  need  to  evaluate  tha  coefflelente 
of  Q44,  q53,  and  2x2x2  Gauea  rule  was  need  to 
Integrate  all  other  term a. 

Since  tha  coefficient  matrix  [K]  depends  on  the 
unkovn  aolutlon  vector  { u) ,  one  should  employ  an 
Iterative  solution  procedure  to  solve  the  finite- 
element  equations*  Bare  wm  use  the  Picard  type 
iterative  technique,  which  begins  with  an  assumed 
displacement  field  (usually,  set  to  taro  to  obtain  the 
linear  aolutlon)  to  compute  [K]  at  the  beginning  of 
the  first  Iteration*  In  subsequent  iterations,  the 
solution  obtained  from  the  previous  Iterations  Is  used 
to  compute  (K|.  The  iteration  Is  continued  until  the 
solutions  obtained  in  two  consecutive  iterations 
differ  by  a  preassigned  error  margin  (say,  10**) • 

Discussion  of  the  Bhmsrlcal  Basalts 


Bare  me  present  results  of  the  linear  analysis  of 
three  test  cases  that  have  been  analysed  by  Pagans 

[20]: 

1*  A  symmetric  three-ply  square  (a  •  b)  laminate 
with  direction  1  of  material  principal  sms 
coincides  with  the  x^-dlrectlon  In  tha  outer 
layers  sod  direction  2  Is  parallel  to  the  x«- 
axla  la  the  center  layer*  The  layers  are  of 
equal  thickness* 

2*  The  same  lamination  geometry  as  la  (1)  above, 
except  that  b  ■  3a* 

3*  A  square  (a  •  b)  sandwich  plate  with  the 
thickness  of  the  face  sheets  equal  to  h/10, 
where  h  is  tha  total  thickness  of  the 
laminate. 


In  all  three  cases,  the  loading  is  assumed  to  be 
sinusoidal  (the  coordinate  system  is  taken  at  the 
center  of  the  plate). 

VX1  WX2 

fc3  2  **xi»x2*  *  «0  °°S  “  c0*  05) 

and  tha  boundary  conditions  ere  of  the  alaply 
supported  type  which  allow  normal  displacement  on  the 
boundary,  but  prevent  tangential  displacement*  For  a 
quartar  plats  tha a a  imply 

At  Xy  -  a/2:  -  t^  •  0 


At  Xj  •  0:  Uj  »  0,  ty  *  0,  t#  •  0 
At  %2  ■  b/2:  Uj  -  u^  •  ty  *  0 

non 

At  Xj  -  0:  “j  "  tx  *  cg  *  0* 


(16) 


Tha  matarial  of  the  laminae  In  cnees  1  end  2  end 
face  sheets  in  case  3  la  assumed  to  have  the  following 
values  for  the  material  coefficients : 


Ej  -  1.72*  *  10®  M/*2  (25  *  10*  p.l) 

S2  »  Ej  •  6.89  x  10*  Ml/*2  (10*  p.l) 

C12  "  c13  “  3,45  *  10*  kH/"2  <°*3  *  10< 
«23  -  13.78  x  10*  Ml/*2  (0.2  x  10*  p.1) 


P*D 

(17) 


The  propertlee  of  the  core  materiel  in  cnee  3  ere 
given  by 

Et  •  Ej  -  0.273  x  10*  Ml/*2  (0.0*  x  10*  p.l) 

8}  -  3. *5  x  10*  Ml/*2  (0.3  x  10*  pal) 

Q|3  «  633  *  0.413  *  10*  Ml/*2  (0.06  x  10*  p.l) 
Cu  -  0.11  x  10*  Ml/*2  (0.016  x  10*  p.1) 


V12  “  W31  "  W32 


0.25 


(18) 


The  following  normalisations  are  used  to  present  the 
results: 


<VV‘l2>  “  — ?  (0l'°2*ff12> 
V 


<®23,013>  "  q^T  <<,23,913> 

B.h2  (E,h3  u,)102 

- - j  u  ,  - ; - 

V  *o* 


(19) 


Table  1  contains  e  comparison  of  the  present 
finite -element  solution  (for  various  meshes)  for  the 
transverse  center  deflection  end  stresses  with  the 
exact  solution  of  Pagsno  [20]  for  Case  l*  As  ehs  mesh 
Is  refined,  the  finite-element  solution  converges 
toward  the  exact  solution*  Further  mesh  refinements 
wars  not  done  due  to  tho  limitations  on  computational 
reoources  (l.a*,  storage,  computational  time,  etc.). 
From  an  examination  of  tha  results  one  can  aas  that 
tha  flnlta-elemeut  results  agree  better  with  the  exact 
solution  for  Increasing  values  of  spun-to-depth  ratio 
(s/h)*  The  transverse  shear  stresses  are  less 
accurately  predicted  than  tha  other  stresses* 


Similar  information  la  prasemtmd  in  Tables  2  end 
3  for  Cases  2  sad  3,  respectively*  Figure  2  contains 


a  plot  of  tha  triMftrN  ihtar  itrasi  distribution 
over  tha  thickness  for  a/h  •  10*  Thn  plot  agraai , 
qualitatively,  vary  nail  with  tha  plot  presented  In 
Figure  7  of  Reference  20  for  a/h  -  4, 
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Table  1.  Compart* on  of  the  finite-* lament  result*  with  the  exact  solution  of  Case  1 
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Figure  1.  Transformation  of  the  master 

element  to  an  arbitrary  element 
of  a  finite-element  mesh 
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Table  2.  Comparison  of  the  finite-element  results  with  the  exact  solution  of  Case  2 
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Table  3.  Comparison  of  the  finite-element  results  with  tha  exact  solution  of  Case  3 
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